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Motivated by the importance of geometric information in real systems, a new model for long- 
range correlated percolation in link-adding networks is proposed with the connecting probability 
decaying with a power-law of the distance on the two-dimensional(2D) plane. By overlapping it 
with Achlioptas process, it serves as a gravity model which can be tuned to facilitate or inhibit 
the network percolation in a generic view, cover a broad range of thresholds. Moreover, it yields 
a set of new scaling relations. In the present work, we develop an approach to determine critical 
points for them by simulating the temporal evolutions of type-I, type-II and type-Ill links(chosen 
from both inter-cluster links, an intra-cluster link compared with an inter-cluster one, and both 
intra-cluster ones, respectively) and corresponding average lengths. Numerical results have revealed 
objective competition between fractions, average lengths of three types of links, verified the balance 
happened at critical points. The variation of decay exponents a or transmission radius R always 
shifts the temporal pace of the evolution, while the steady average lengths and the fractions of 
links always keep unchanged just as the values in Achlioptas process. Strategy with maximum 
gravity can keep steady average length, while that with minimum one can surpass it. Without the 
confinement of transmission range, I — > oo in thermodynamic limit, while I does not when with it. 
However, both mechanisms support critical points. In two-dimensional free space, the relevance of 
correlated percolation in link-adding process is verified by validation of new scaling relations with 
various exponent a, which violates the scaling law of Weinrib's. 
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I. INTRODUCTION 



Correlated percolation jl|-|14{ is a useful theoretic model 
in statistical physics. It provides us with fundamental un- 
derstanding of spread processes of message, disease and 
matter in nature and society. Linking probability be- 
tween any two nodes in it takes the form of p(r) ~ r~ a , 
where r is d-dimcnsional distance between the nodes, and 
a is a positive real number, namely, distance-decay ex- 
ponent of links. Weinrib and Halperin[lj| analytically 
studied whether the correlations change the percolation 
behavior or not. Weinrib pointed out, for a < d, 
the correlations are relevant if av — 2 < 0, where v is the 
percolation-length exponent for uncorrelatcd percolation; 
while for a > d the correlations arc relevant if dv — 2 < 0. 



It is a generalization of the Harris criterion[17| appears 
earlier. Recently, network models referring to correlated 
percolation have gradually appeared (Tl|. 

Achlioptas process (AP) [18| for link- adding networks, 
which is an attractive topic at present (l9H38j . could be 
viewed as a new kind of correlated percolation if we 
put all nodes uniformly on a two-dimcnsional(2D) plane. 
Starting from a set of isolated nodes, two candidate links 
are put to nodes randomly at every time step, but only 
the link with smaller product m^m^- is retained, where 
mi (or rrij) is the mass (the number of nodes) of the clus- 
ter that node i(or j) belongs to, which is called Product 
Rule (PR). A link chosen with PR from both inter-cluster 
candidates is called a type-I link. When two candidate- 
links are of different types, i.e., one is an inter-clustcr 
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link, the other is a intra-cluster one, always the later is 
retained, and it is called a type-II link. While for both 
intra-cluster ones, the retained link is arbitrarily chosen 
no matter they are in the same or different clusters, and 
it is called a type-Ill link. Generally speaking, in the 
way of AP, network percolation is inhibited, which post- 
pones the appearance of the threshold T c at which a gi- 
ant component G starts to grow, and results in a sharp 
growth of G called an explosive percolation. In our point 
of view, if we put AP on a 2D plane, it gives rise to a 
new mechanism of long-range correlation for the nodes 
based on co-evolutionary growing masses of components 
they connected. The selective rule for topological links 
relies on mass-product instead of 2D geometric length of 
them, which prevents the property exhibited in the pre- 
vious correlated percolation. And correlation feature in 
AP-type of percolation has not been revealed up till now. 

A recent model of ours[3Ej based on the observation 
of phenomena in different real systems [40l - [45| describes 
another kind of correlated percolation in growing net- 
works, which could be viewed as a overlapping of tra- 
ditional correlated percolation with AP in a 2D space. 
The link-occupation function in the model takes the form 
p(r) ~ miirij/rfj which looks like Newton's gravity rule. 
It resumes the classical Erdos-Rcnyi(ER) random graph 
model when exponent a — ¥ oo, and it gives another ex- 
treme of AP when a — > 0. Different properties of such 
kind of new correlated percolation are expected, since ER 
random graph grows without any bias, AP takes a strong 
bias to inhibit network percolation independent of geo- 
metric distance, while the new model with gravity-like 
rule have some distance-related relax on such bias, which 
produces a new type of correlated percolation. 

In this paper, we report the simulation results on ob- 
jective competition between type-I, II and III links in 
both gravity model[43|, H|| and AP model, and we point 
out a new mechanism to support critical points, which 
bears the scaling relations revealed in our recent work. 
Different saturation effect is manifested, which distin- 
guishes it from the traditional correlated percolations. 



II. MODEL 

Suppose N isolated nodes are uniformly scattered on 
a 2D plane. For convenience of calculating distance, the 
plane is discretized with a triangular lattice, each mini- 
mal edge with the length of two units for the convenience 
of algorithm. Each vertex of the triangles is occupied by 
a node so that we exclude all possible biases except link- 
adding rules. For any two pairs of nodes i and j possibly 
with the same product m^m.,-, a type-I link connects the 
the pair with longer distance if both the links' ends hit 
the nodes belonging to different clusters; while a type- 
II link connects the nodes inside the same cluster if the 
other one is an inter-cluster link; a type-Ill one connects 
arbitrarily chosen pair of nodes if both candidate links 
arc intra-cluster ones. 



Parallel to PR, we pick randomly two pairs and 
(k, I)] of nodes in the plane at every time step. For the 
pair (and for (k,l) likewise), we calculate the gen- 

eralized gravity defined by gij = miTrij/rfj, where 
is the geometric distance between i and j, and a is an 
adjustable decay exponent. Once we have gij and gki, 
we have two choices in selecting which pair to connect. 
For the case of the maximum gravity strategy (we call it 
G max ) we connect the pair with the larger value of the 
gravity, e.g., the link is made if <?y > gu and the 
link (k,l) otherwise. We also use the minimum gravity 
strategy (G m i n ) in which we favor the smaller gravity pair 
to make connection. The two strategies, G max and G m i n , 
lead the link-adding networks to evolve along the oppo- 
site paths of percolation processes. Generally speaking, 
G max facilitates the percolation process, whereas G m - m 
inhibits it. All such generalized gravity values are cal- 
culated inside the circular transmission range with the 
radius R centered at one of nodes i and j a s the speaking 
node [4(| in a mobile ad hoc network [4li l42l |45| . For 
the different limits of parameters R and d, we have three 
cases in the model. Case I: With the transmission range 
R — > oo, we have a generalized gravitation rule which is 
an extension [45[ of widely used gravitation model [43| 
(a = 1) with the tunable decaying exponent a. Case II: 
With the exponent a = 0, we assume that node pairs 
can be linked with PR topologically inside the transmis- 
sion range with a limited radius R. Case III: With both 
limited values of radius R and exponent a, we have the 
gravity rule inside the transmission range. It can describe 
the communication or traffics with constrained power or 
resources. 

For case I and case III in the model, three scaling rela- 
tions have been found with large scale simulations. When 
strategy G max is adopted in 2D free space(case I), we 
have 

C-a^fita^) (1) 

where t = (T — Tq)/Tq is dimensionless time-step with 
To = 0.78, a is the decay exponent of connection proba- 
bility, a = 0.01, e = 0.20, and f(x) is a universal func- 
tion. When strategy G m i n is adopted inside the trans- 
mission range with radius i?(case III), we have 

C ~ (a/a )- 9 h[t{a/a f) (2) 

for certain parameter ranges of a and R, where t = (T — 
T ) /T , T = 1.0, 9 = 0.005, (j) = -0.50, a = 0.5, and 
h{x) is a universal function. In addition, when strategy 
Gmax for case III is adopted inside transmission range 
defined by R, we have another scaling relation 

C ~ R- s H{tff>) (3) 

for R > 3, where p= (R- R )/R , R = 2, 77 = -0.10, 
5 = —0.005, To = 1.0 and H(x) is a universal function. 

To understand three scaling relations above, we should 
look into the mechanism of the evolution processes under- 
lying C{T). To sec what happens in such critical points 
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FIG. 1: (color online)Evolution of fractions of type- 1, type-II 
and type-Ill links in Achlioptas process on 2D plane. 



FIG. 2: (color online)Evolution of average lengths of type-I, 
type-II and type-Ill links in Achlioptas process on 2D plane. 



To, and what are particular of them in certain link-adding 
processes, we count the temporal link fractions F(T), and 
calculate the average lengths I of links which is defined as 
the summation of all lengths of links for a certain type 
over its number in a window AT = 20 time-steps. By 
observation of the time-dependent behaviors of fractions 
of type-I, II and III links, new properties were found out 
for our gravity-like model together with AP producing 
explosive percolations. 



III. 3. SIMULATION RESULTS 

All simulations are carried out on the triangular lattice 
of the size N = L x L with L = 32, 64, 128 and 256, 
respectively. We simulate either of strategy G max or G m ; n 
for either case I or III. The total number of links equating 
to that of time-steps is divided by N, which is defined 
as T. The mass of the largest component divided by 
N makes up the observable C, the node fraction of the 
largest component. All results presented in this work 
are obtained from 5000 different realizations of network 
configurations with L = 128 if not specially indicated. 

Inspired by Cho and Kahng's [33j work and a referee 
of ref.(39j. we have gone further by calculating fractions 
of three types of links and arithmetic average lengths of 
their links. Our attention was pointed at AP first. In 
Fig.l we illustrated the evolution of fractions of 3 types 
of links. Just at the threshold the fraction of type-l(Fj) 
links has a sharp drop-down, meanwhile that of type-II 
(F H ) shoots up, crossing Fj at T c = 0.888. A little af- 
ter it, Fj crosses with growing fraction of type-Ill links 
{Fin) at the level Fi = Fm = 0.25, while Fn gets its 
summit(T/7 = 0.5) at the same point, which has not been 
concerned by previous works. However, it is this property 
that pervades all cases in the present correlated percola- 



1.0 

0.8 

0.6 

0.4 

0.2 

0.0 






















a 










0.2 1 










— 0.2 II 










0.2 III 










— 0.5 1 - 










0.5 II 




\ 














— 0.5 lll^_ 
1.0 1 ^ 


















Igg 1.0 II 










1.0 III 










2.0 I 










gjj 2.0 II I - 




ft ^"^^ 






— 2.o iiiS 

3.0 I 


















— 3.0 II — 




+fn -f- 










11' S 






- 3.0 III 










^-T 5.0 I 

5.0 II — 










: 5.0 III 











0.5 



1.0 



1.5 



2.0 



2.5 



FIG. 3: (color online)Evolution of fractions of type-I, type-II 
and type-Ill links with G max and probability decay exponent 
a = 0.2, 0.5, 1.0, 2.0, 3.0ond5.0 in case I (R -> oo). 



tion. In Fig. 2, the average lengths of type-H(Zjj) merges 
that of type-III(Z///) at T c after an abrupt growth, In 
starts to grow earlier than //// . The level of I for both of 
them keep invariant for T > T c , while li starts to decrease 
from T c . We see from both the figures that in explosive 
percolation the system undergoes a sharp transition from 
a type-I link dominant phase into a type-II and III dom- 
inant phase at T c [33[. Besides, average lengths undergo 
a parallel transition at the same point. Actually, 3 levels 
of I go to infinity in dynamic limit from finite size scaling 
transformation(not shown). 

Now we turn to the fraction of 3 types of links in case I 
of the present gravity model [3{|. With strategy G max in 
free 2D space, we have scaling relation (1) for distance- 
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FIG. 4: (color online) Evolution of average lengths / of type- 
I, type-II and type-Ill links with Gmax and probability decay 
exponent a = 0.2, 0.5, 1.0, 2.0, 3.0and5.0 in case I. 
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FIG. 5: (color online) Node fraction C(T) of the largest com- 
ponent with strategy G max and the same parameters in Fig. 3. 



decay exponent a € [0.2,2.0]. Fig. 3 shows the evolution 
of fractions corresponding to it. Curves for Fi cross those 
of F n with a = 0.5, 1.0 and 2.0 around T = 0.78 quite 
clearly, even those with a = 0.2 barely cross near it. 
But fractions for a = 3.0 and 5.0 shift the cross point 
rightward obviously. On the other hand, the curves for 
Fj cross Fin at T = 1.0 for almost all values a except 
a = 0.2. To determine which one would be the candidate 
of another critical point, we cast ourselves on the assis- 
tance from the observation of average lengths of links. In 
Fig.4 T U {T) merges l n (T) for a = 3.0 and 5.0 at T = 1.0, 
separating themselves from others. However, almost all 
other ones collect at T = 0.78, which gives hint to us 
for Tq. (Here a better resolution is needed in further cal- 
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FIG. 6: (color online) (a)Susceptibility x(T) of the system 
with Gmax in Case I. L = 32,64, 128 and 256.(b)Percolation 
thresholds T c for a = 0.2,0.5,0.8,1.2,2.0,3.0 and 5.0 with 
Gmax in Case I. 



culation). Correspondingly, in Fig. 5 for scaling relation 
(1) with exponents e = 0.2, a = 0.01, and Tq = 0.78, 
C(T) for a G [0.5,2.0] collapse into the universal func- 
tion very well, with that for a = 0.2 barely collapsing 
onto it. But those C(T) for a = 3.0 and 5.0 do not be- 
have well in collapse. The separation from others at the 
turning middle part indicates the deviation of their Tq 
from 0.78 with which others share. In the description 
of the average lengths for case I with G m ax, simulated 
results T(T) in Fig.4 with all values a demonstrate the 
same steady level {I ~ 131.50 for L = 128). Variation 
of parameter a only shifts starting points of up-growing 
Fjj and Fjjj as a increases. The saturation effect of 
large decay exponents (a = 3.0, 5.0) appears clearly and 
is shown by dash lines, which demonstrates the inheri- 
tance from traditional correlated percolationflll]. In this 
case with G m ax, special level of fractions at cross point 
of Fj{T) and F ln (T) keeps 0.25 just as in AP without 
any distance-decay included, so does Fn at hiking its 
summit. 

As in a usual way, we determine the critical point T c of 
percolation by observation of tips of susceptibility x(^) 
(Fig. 6). Comparing T c in Fig. 6 with l(T) in Fig.4, we 
find that these T c approximately hit the horizontal co- 
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FIG. 7: (color online) Size-dependent average lengths I of 
type-I, type-II and type-Ill links with strategy G max in case 
I. Parameters are the same as in Fig. 6 



ordinates of middle point of growing fraction of type-II 
links, which means that T c is the transition point from 
the inter-clustcr-link dominant phase to the intra-cluster- 
link dominant phase. Besides T c , we have another (sub)- 
critical point To which is in certain range independent 
of decay exponent a in gravity model, and To indicates 
the balance between the fractions of type-I and type-II 
links, yielding a new scaling behavior of C(T, a) in for- 
mula (1) not revealed by previous works. Moreover, the 
steady level I = 131.5 is always a— independent, type- 
independent which takes the inherited value of that in 
AP. Actually, 131.5 is the value for L=128 only. We have 
size effect since a free boundary condition instead of a 
periodic one is adopted. The finite size effect is shown in 
Fig. 7 which gives that I ~ L , i.e., 



T-N 1 ' 2 



(4) 



where N is the number of nodes on the 2D plane. 
Hopefully it goes towards infinity in thermodynamic 
limit. However, the finite size effect of C(T,N) (Fig. 8 
for an example) is not strong enough for us to identify 
the scaling exponents v, f3 as usual. Therefore, we can 
check the validation of scaling laws presented by Weinrib 
[l(| for correlated percolation in the present model only 
by rescaling susceptibility x(o, T). Fig. 9 illustrates the 
results of it for examples a = 0.5 and 2.0, respectively. 
We have scaling relation 



X ~ N l/v G{{T - T c )N 1/v ) 



(5) 



where \jv = 0.3, j/v = 0.75, and G{x) is a universal 
function. With these values of scaling exponents, the 
scaling law vi ong = 2/ a is not applicable to the present 



1.0 
0.8 
0.6 
0.4 
0.2 



0.0 - 



■ 

- 




a = 0.5 






N 




32 


/ 


64 


; J 


- 128 




256 


i i i i i i i i i l i i i i i i i 





0.0 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



0.5 1.0 
T 



1.5 







a = 2.0 






N 




32 




64 




128 




256 






' ■ 





0.0 0.5 



1.0 

T 



1.5 



2.0 



FIG. 8: (color online)Node fraction C(T) of the largest com- 
ponent with strategy G max in Case I. (a)a = 0.5; (b)a = 2.0. 
N = 32, 64, 128 and 256. 



model, where vi ong is correlation-length exponent in the 
long-range case. 

With power-law form for the correlation function g(r), 
Weinrib had derived the extended Harris criterion: the 
long-range nature of the correlations is relevant if av— 2 < 
0, which means the correlations change the percolation 
critical behavior. It has been violated since now they all 
behave differently from traditional short range percola- 
tion in a 2D triangular lattice (y = 4/3) and the cor- 
relations are relevant no matter av — 2 is less(Fig.9(a)) 
or larger(Fig.9b) than zero. This is because in strat- 
egy G max we have overlapped the power-law correlation 
function g(r) with AP which is another kind of autocor- 
relation process with positive feed back effect of mass- 
growing. 

For the strategy G m i n which prefers smaller gravity, it 
tends to retain a longer link under the comparison of the 
same product of masses rmmj. That is to say, long range 
links have predominance. Evolution of Fi(T) and Fu(T) 
links do not cross at any common point. In Fig. 10, the 
cross points for Fj and Fjj shift leftward from T c of AP 
as decay exponent a increases, which means that G m i n as 
a correlated percolation mechanism weaken the explosive 
effect caused by AP. But the starting position of fraction- 
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FIG. 9: (color online) Scaling of susceptibility x(T) in case I 
for (a) a = 0.5; (b) a = 2.0. N = 32, 64, 128 and 256. 
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FIG. 10: (color online)Evolution of fractions of type-I, type-II 
and type-Ill links with Gmin and probability decay exponent 
a = 0.2, 1.2, 2.0 and 5.0 in case I. L = 256. 




FIG. 11: (color online) Evolution of average lengths of type-I, 
type-II and type-Ill links with Gmin and probability decay 
exponent a = 0.2, 1.2, 2.0 and 5.0 in case I. L = 256. 



II still provides hints of thresholds T c . However, it is 
hard to locate a common cross point for .Fj and Fm in 
a range of a. Therefore, we have no scaling relation for 
them. The steady average lengths of links (Fig. 11) have 
the same level of G max and AP cases, and size-effect (not 
shown) also tells the divergence of I, but all of them do 
not tell any possible hint for critical points. 

Generally speaking, strategy G m i n facilitates longer 
links for certain geometric distribution of clusters or 
nodes. In the evolution, strategy G m i n emphasis the as- 
signments for different types of links, encourage longer 
and intra-cluster links. Humps above the steady Fjj 
and Fm implies out-of-pace growing of link lengths of 
Fj. That is, Fj surpasses the growing speed of the giant 
component. Here, it is the geometric distance-dependent 
strategy that makes G m i n alleviate effect of AP. With 
smaller a (e.g., a=0.2 in Fig.ll)the strategy has the op- 
portunity to exhaust long links before T c ; while with mid- 
dle values of a (e.g., a = 1.2 and 2.0) it may take longer 
time to exhaust them. However, with too large a (e.g., 
a = 5.0), G m i n fall off quicker than the natural dimen- 
sion, we can only see the AP-type short-range effect of 
saturation. In this limit, i.e., a < 3.0, percolations are no 
longer relevant, which causes saturation of curves C(T) 
in Fig. lb [13 • However, we should not expect the short- 
range percolation exponent v = 4/3 of correlation length 
here for 2D triangular lattice, since AP has been included 
in the strategy G m i n . 



The distinct feature in case III for both G„ 



and 



Gmax is that the candidate links are selected not only by 
comparing gravities, but also constrained inside a trans- 
mission range R (r = 2R in geometric distance), which 
ruins the effect comes from the divergence of average link 
lengths. It is well known that all possible singularities at 
critical points come from the singularity of correlation 
length. However, here no length could goes to infinity in 



7 















3 














_ U.i! _ 
























R=4 




2 III 












0.5 1 














0.5 11 


















: vj'. \\\\ 










» 1.0 1 f- 




\\ \\\\\ ■ 










•! 1.0 II 

1.0 III 

1 .5 I - 














1.5 11 














1.5 III 




1 








2.0 I 












2.0 II 












2.0 III 












3.0 I 




/ It IS. 








2 3.0 II 

g 3.0 














5.0 I 












5.0 II 












5.0 III 

. 10.0 I 












■ -10.0 II - 












10.0 III 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



T 



FIG. 12: (color online)Evolution of fractions of type-I, type-II 
and type-Ill links with strategy Gmin and probability decay 
exponent a = 0.2, 0.5, 1.0, 1.5, 2.0, 3.0, 5.0 and 10.0 in case III. 
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FIG. 14: (color online) Scaling of node fractions C(T) with 
strategy Gmin, a = 1.0, 2.0 and 3.0in Case III. R — 4. 
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FIG. 13: (color online) Evolution of average lengths of type-I, 
type-II and type-Ill links with strategy Gmin and probability 
decay exponent a = 0.2, 0.5, 1.0, 1.5, 2.0, 3.0, 5.0 and 10.0 in 
case III. R = 4. 



any way, which ruins possible common cross point relies 
on the balance between Fj and Fjj as in the 2D free space 
(case I of the present model) and induces the possibility 
to yield novel scaling relation other than any previous 
ones. For possible critical point, we seek help from the 
evolution of link fractions of 3 types. Fig. 12 shows the 
behaviors of Fj, Fjj and Fm with G m i n for all simulated 
distance-decay exponents a and R — 4. 

Correspondingly, Fig. 13 shows the behaviors of ij, In 
and Ijjj the the same set of parameters. Critical point 
T = 1.0 (Tq = 0.99 to be precise) distinguishes itself 
from others by intuitive observation. The cross point for 
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FIG. 15: (color online)Pseudo-scaling of node fractions C(T) 
with strategy Gmin, a = 0.2, 0.5, 1.0, 1.5, 2.0, 3.0, 5.0 and 10.0 
in Case III. R = 4. 



a = 5.0 and a = 10.0 go rightward from which others 
share, which means they could not share the same To 
hence the same scaling relation with other decay expo- 
nents. Besides, in rescaling process for C(T, a), curves for 
a = 0.2 and a = 0.5 failed in collapse, because smaller 
transmission range R inhibits the effect of slower (long- 
range) decay for connection probability. The rescaled 
function C{t) for R = 4 is shown in Fig. 14. It seems that 
we could go further with <p = —1.5 to include more expo- 
nents a in the scaling as shown in Fig. 15. However, it is 
meaningless in physics due to above mentioned reasons. 
Actually, scaling behaviors are R-depcndent, but expo- 
nents (j) and 9 need not to vary. The variation of R only 
shifts Tq as illustrated in Fig. 16 for R = 8. Therefore, we 



8 




5.0 1 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

T 



FIG. 16: (color online)Evolution of fractions of type-I, type-II 
and type-Ill links with strategy G m in and probability decay 
exponent a — 0.2, 0.5, 1.0, 1.5, 2.0, 3.0 and 5.0 in case III. R = 
8. 
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FIG. 18: (color online)Evolution of fractions of type-I, type-II 
and type-Ill links with strategy Gmax and probability decay 
exponent a — 5.0 for R = 3, 4, 6 and 8 in case III. 
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FIG. 17: (color online) Evolution of fractions of type-I, type-II 
and type-Ill links with strategy Gmax and probability decay 
exponent a = 2.0 for R — 4,6 and 8 in case III. L = 256. 500 
realizations of network configurations. 



FIG. 19: (color online) Evolution of average lengths of type-I, 
type-II and type-Ill links with strategy Gmax and probability 
decay exponent a = 2.0 for R = 4, 6 and 8 in case III. L = 256. 
500 realizations of network configurations. 



keep (j) = —0.5, 9 = 0.005 for all values of R with G m j n , 
but take T = 0.99 for R = 4, T = 0.92 for R = 8, 
and so on. The humps above the steady level of In and 
In i (all independent of exponents a) come from similar 
mechanism as in free 2D space but now at much lower 
level constrained by transmission radius R, and they are 
independent of size L of the system. 

Scaling relation (3) for case III with strategy G max 
inside transmission range with radius R is checked for 
various decay exponents a and for different sizes (L=32, 
64, 128 and 256). Its validation is independent of size 
L simulated. In Fig. 17 and Fig. 18, the evolution of Fj, 



Fji and Fjjj for both a = 2.0 and a = 5.0 behave much 
similarly. Fj and Fjj cross at the level a little bit lower 
than 0.45, while Fi crosses Fm at the level 0.25, which 
keeps the same as in all previous cases. The changes of 
exponent a and R only shift fractions along horizontal 
direction of figures, i.e., to change starting points and 
growing/dropping speed instead of levels of them. How- 
ever, To = 1.0 keeps as their common fixed point for Fj 
and Fm to cross. In Fig. 19 and Fig. 20, Fi and Fm for 
both a = 2.0 and a = 5.0 arrive at the same level hitting 
To = 1.0, which distinguishes this point from totally 3 
cross points, and makes up a candidate of critical point 
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FIG. 20: (color online) Evolution of average lengths of type-I, 
type-II and type-Ill links with strategy Gmax and probability 
decay exponent a — 5.0 for R — 3, 4, 6 and 8 in case III. 
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FIG. 21: (color online) Scaling of node fractions C(T) with 
strategy Gmax, a = 2.0 for R = 4, 6 and 8 in Case III.L = 256. 
500 realizations of network configurations. 



T for scaling relations. The scaling exponents r\ = —0.10 
and 6 = —0.005 have been checked for lower values of pa- 
rameter a(0.5 < a < 3.0, Fig. 21). However, for a = 5.0, 
we have to choose a new set of exponents: r\ = —0.25 and 
S = —0.01 (Fig. 22). It is not strange that steady levels of 
I keep unchanged for certain R, independent of a or L, 
just as that with G m i n in case III. I inside a circle defined 
by R can not go to infinity under any circumstance, but 
still support a critical point, which distinguishes corre- 
lated percolation in case III from case I and traditional 
models. It deserves further investigation. 
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FIG. 22: (color online) Scaling of node fractions C(T) with 
strategy Gmax, a = 5.0 for R = 3, 4, 6 and 8 in Case III. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have proposed a new network model 
of correlated percolation in which geometric distance- 
dependent power-law decay connection probability over- 
laps Achlioptas process to form a gravity model. It can 
be tuned to facilitate or inhibit percolation with strategy 
Gmax or Gmin, cover a wide range of thresholds T c , yield 
a set of new scaling relations. And it provides a scheme 
for better description of practical processes in complex 
systems. 

We have developed a new approach to find out can- 
didate critical points with physical meanings other than 
that of traditional ones. There are objective competition 
and balance between type-I and type-II, type-I and type- 
Ill links, meanwhile, competition of average lengths be- 
tween type-II and type-Ill links. Along this line threshold 
T c is found to overlap the balance point between factions 
Fi and Fjj in the explosive percolation of Achlioptas 
process, and the steady average lengths of three types 
of links are all divergent to infinity in thermodynamic 
limit. The percolation is indeed a transition from type- 
I link dominant phase to type-II and type-Ill dominant 
phase. 

By observing evolutions of fractions of type-I, type-II 
and type-Ill links, a candidate critical point can be cho- 
sen combined with the message on evolutions of average 
lengths of them. With strategy G max in 2D triangular 
lattice, fraction Fj get balance with Fn , makes up a crit- 
ical point Tq which supports scaling relation (1) in case 
I of the model. With strategy G m in inside certain trans- 
mission range with radius R, a duet balance exists for Fj 
and Fiji meanwhile In and lm, makes up another crit- 
ical point To which supports scaling relation (2) in case 
III. With strategy G m in and certain range 0.5 < a < 3.0 
of decay exponent a, again a duet balance exists for Fj 



10 



and Fin meanwhile In and lm, makes up another crit- 
ical point T which supports scaling relation (3) for a 
mini-scale of R in case III. This approach serves an as- 
sistant tool in seeking critical points of order parameter 
C(T) which is usually not easy to determine in an intu- 
itive way. 

In numerical calculations, besides percolation thresh- 
old T c , two fixed points, To = 0.78 and To = 1.0 emerge as 
distinct points not only for special temporal crux but also 
for unchanged levels of Fj, Fjj and Fjjj inherited from 
AP, which is expected to be further proved in analytical 
ways. However, they have different physical meanings. 
The former corresponds to a divergent average length 
of links, while the later corresponds to confined average 
lengths by transmission range R, which distinguishes it- 
self from traditional critical points in percolations. 

Correlated percolations are relevant since long-range 
correlation drastically changes the critical properties. 
The validation ranges of decay exponents a with vari- 
ous strategies in different cases define the relevance of 
correlation. They have demonstrated novel scaling rela- 
tions different from traditional 2D short-range percola- 
tion in triangular lattice. The intervention of distance- 
dependent power-law decay ingredients alleviates the ex- 
plosive effect of percolation transition by horizontal ad- 
justment of evolutions along temporal axis, separates Tq 



from T c , while the overlapped AP included in the present 
gravity model always conquers the vertical levels of three 
fractions and average lengths, which are found neither in 
traditional correlated percolations of continuities in 2D 
space nor in complex networks. Moreover, the node frac- 
tion C(T) of the largest component, fractions F(T), and 
average lengths T(T) of three types of links all show sat- 
uration phenomena as pointed out by Weinrib but with 
different values of exponent of a since now AP overlaps in 
the present gravity model. And scaling law of Weinrib is 
no longer obeyed according to the evidence of numerical 
results of average-length exponents. 
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